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Dispersive shock waves dominate wave-breaking phenomena in Hamiltonian systems. In the 
absence of loss, these highly irregular and disordered waves are potentially reversible. However, no 
experimental evidence has been given about the possibility of inverting the dynamics of a dispersive 
shock wave and turn it into a regular wave-front. Nevertheless, the opposite scenario, i.e., a smooth 
wave generating turbulent dynamics is well studied and observed in experiments. 

Here we introduce a new theoretical formulation for the dynamics in a highly nonlocal and defo- 
cusing medium described by the nonlinear Schroedinger equation. Our theory unveils a mechanism 
that enhances the degree of irreversibility. This mechanism explains why a dispersive shock cannot 
be reversed in evolution even for an arbitrarirly small amount of loss. 

Our theory is based on the concept of nonlinear Gamow vectors, i.e., power dependent gener¬ 
alizations of the counter-intuitive and hereto elusive exponentially decaying states in Hamiltonian 
systems. We theoretically show that nonlinear Gamow vectors play a fundamental role in nonlinear 
Schroedinger models: they may be used as a generalized basis for describing the dynamics of the 
shock waves, and affect the degree of irreversibility of wave-breaking phenomena. Gamow vectors 
allow to analytically calculate the amount of breaking of time-reversal with a quantitative agreement 
with numerical solutions. 

We also show that a nonlocal nonlinear optical medium may act as a simulator for the experimental 
investigation of quantum irreversible models, as the reversed harmonic oscillator. 


I. INTRODUCTION 

One of the fundamental problems in physics is the de¬ 
scription of irreversible processes. Irreversible processes 
are denoted by exponentially decaying observables asso¬ 
ciated with a preferential arrow of time. However, ex¬ 
ponential dynamics of the wavefnnction can be excluded 
from first principles in leading theoretical models, such as 
the Hilbert space (HS) formulation of quantum mechan¬ 
ics (QM). In this formulation, there are theorems stating 
that the solution of the Schrodinger equation cannot have 
an exponential evolution. [I, 2] 

To explain damping, one introduces loss in a phe- 
nomenogical way, accounting for the coupling with the 
environment. [3, 4] However, exponential time decay of 
observable quantities, which signals irreversible dynam¬ 
ics, is routinely detected in experiments, even when the 
role of the environment is negligible. 

This is known as the probability problem and is common 
to quantum and classical systems. [5] 

A paradigmatic example of the issues related to irre¬ 
versibility in classical systems is found in nonlinear wave 
propagation. In regimes allowing for shock waves, singu¬ 
lar solutions may be excited by smooth initial conditions 
that generate an evolution dominated by non-regular, 
highly disordered, waves (wave-breaking). [6] The evo¬ 
lution of smooth initial conditions towards a highly dis¬ 
ordered and incoherent regime is intuitive. However, the 
reverse process, albeit possible, is strongly puzzling, and 
against well accepted ideas about entropy. 

Dispersive shock waves (DSW) and wave-breaking 
in Hamiltonian models therefore question reversibility. 
They seem to indicate a preferential arrow of time that 


should be excluded from first principles, because no cou¬ 
pling with the environment is foreseen. [7-13] 

Here we pose the following question: is there a con¬ 
nection between shock waves and the problem of the ir¬ 
reversibility of the wavefnnction? These two problems 
are seemingly unrelated, but if one carefully looks at the 
dynamics of a shock wave (as detailed below), one realizes 
that exponential dynamics can be identified. We show in 
the following that the specific case of a nonlocal nonlin¬ 
earity enables to unveil the link and lay the foundations 
for further developments. 

The starting point of the analysis is that one can cir¬ 
cumvent the absence of exponentially decaying solution 
for the wavefnnction by enlarging the Hilbert space for in¬ 
cluding non-normalizable wavefunctions that decay expo¬ 
nentially in evolution. This is known as the rigged Hilbert 
space (RHS) formulation of quantum [14] and classical 
mechanics [15, 16]. Exponential dynamics are explained 
in terms of the excitation of states with complex en¬ 
ergy eigenvalues. These states, known as Gamow vec¬ 
tors (GVs), were originally introduced in nuclear physics, 
and model irreversibility in Hamiltonian systems. [17-20] 
Even if GV are not-normalizable, they can be used as 
a generalized basis, as it happens, e.g., for plane waves. 
However, the watermark of GV is a quantized decay rate. 
This quantization is radically different from the standard 
quantization of real eigenvalues, and involve the imagi¬ 
nary part of the GV eigenvalues. Notably enough, so 
far, no experimental evidence of the quantization of de¬ 
cay rates has been reported. We also remark that this 
has nothing to do with quantum many-body theory, or 
out-of-equilibrium systems (see, e.g., [4, 21]). 

In this manuscript, we show that GV can be intro- 
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duced to describe highly nonlinear regimes in nonlocal 
media. We show that the solution of the relevant nonlin¬ 
ear Schroedinger equation in the presence of shock waves 
can be given in terms of GV with power-dependent, quan¬ 
tized decay rates. GV are shown to quantitatively de¬ 
scribe shock waves, in terms of the amplitude and the 
phase of the complex wavefunction and they explain the 
observed exponential dynamics. 

We also demonstrate that the excitation of GV radi¬ 
cally changes the degree of irreversibility of the dynamics. 
Specifically, the theoretical tools arising from the use of 
GV enable to explicitily consider and calculate the de¬ 
gree of irreversibility, defined as the difference between 
the forward and the backward propagated wavefunction. 
This calculation is found to be in quantitative agreement 
with numerical solutions, and shows that even an arbi¬ 
trary small amount of loss produces a breaking of time- 
reversal enhanced by shock waves. 

This manuscript is organized as follows: in section II 
we introduce the leading model; in section III we intro¬ 
duce GVs and the resulting expansion; the specific case 
of the reversed harmonic oscillator is analyzed in sec¬ 
tion IV; section V describes the link between the ground 
state GV, shock waves and nonlocal solitons; section VI 
describes the exponentially decaying dynamics; section 
VII deals with the definition of the degree of irreversibil¬ 
ity and the way it can be calculated by the GV expan¬ 
sion, quantitative comparison with numerical simulations 
is also reported; conclusions are drawn in VIII. 


II. NONLOCAL NONLINEARITY 


We start from the paraxial wave equation for the prop¬ 
agation of an optical beam with amplitude A and wave¬ 
length A in a medium with refractive index no, and linear 
loss length Lioss [R- = (^, ^)] [10] 


2^fc 
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no 


A = -i- 


-'loss 
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A is normalized such that I = \A\‘^ is the intensity, 
Pmks = / IdYi is the power, k = 27rno/A is the 
wavenumber. In (1) An is the nonlinear nonlocal per¬ 
turbation to the refractive index. 


An[|Ap](R) = n 2 y G 2 (R - R')/(R')^I^' (2) 

with G 2 the kernel function, normalized such that 
J G 2 dR = 1. The local Kerr effect corresponds to G 2 = 
(5(R-R'). For an exponential nonlocality G 2 (V, V) = 
G(V)G(V), being G(V) = exp(-|V|/L,zoe)/(2T,/oc), 
with Lnioc the nonlocality length. 

Since the problem is separable in the highly nonlo¬ 
cal approximation (HNA) considered below, we consider 
the case of one transverse dimension X. For a defocus- 
ing nonlinearity (722 < 0), we write Eq.(l) in terms of 
dimensionless variables, letting Wq be the beam waist. 



FIG. 1. (Color online) (a) GV |f“(x)|^ for increasing even 
order n; (b) dx arg[f“(x)] for increasing even order; (c) weight 
of the GV expansion of a Gaussian wavepacket. 




FIG. 2 . (a) Numerical solution of Eq.(3) with P = 10^ and 
cr^ = 10; (b) projection on GV for increasing n for a — 0.3 
and 7 = 8 , continuous lines are from Eq.(3), dots are from 
Eq.(20); (c) as in (b) for 7 = 24. 


X = V/kFo, and 2 : = ZjZ^^ with Z^ = kW^ the diffrac¬ 
tion length: 

where a = Z^/Lioss, '*/’ = AWo/y'PMKS, and {tpltp) = 1, 
brackets denoting the HS scalar product. The asterisk 
in (3) is a convolution, and P = PmksIPref with 
Pref = AV4^2^o|n2|, and K{x) = W^G{xW^) = 
exp(-|x|/cr)/2cr with a = Lnioc/Wo. 


III. GAMOW VECTOR EXPANSION 

Starting from Eq.(3), in the HNA, one obtains an ef¬ 
fectively linear Hamiltonian system: we have [22, 23] 

K * |'0(x)p = r{x), (4) 

and 




( 5 ) 
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with the Hamiltonian 

H=\f + V{x), (6) 

being V{x) = Pk{x) and p = —idx- 

GVs are numerable generalized eigenvectors of H with 
complex eigenvalues with n = 0 , 1 , 2 ,...; they exist 

outside the HS of Lebesgue integrable functions. In the 
RHS, they furnish a generalized basis for a normalizable 
wavepacket, being H(l)^ = 

oo 

Ipix, 0) = y] (/)^ {x) + <f)b{x) = <jp{x) + <j)b{x). (7) 

n=0 

In (7) '0(x,O) is the initial state, the discrete part is not 
the usual spectrum of normalizable bound states with 
real valued eigenvalues (which is not present here), but 
an ensemble of vectors with complex eigenvalues 

En = E^- iTn/2, (8) 

with real part and immaginary part correspond¬ 
ing to the decay rates. [14, 18, 24] 

05 is the background component that is due to the part of 
the continuum spectrum of H not containing GVs. The 
key difference between 0 ^ and 05 lies in the fact that, 
upon evolution, the former decays exponentially, while 
exponential decay is not allowed for 05. We have 

z) = e~^^^(tP{x) = Y,n 

During evolution GV cannot emerge from the back¬ 
ground 05 , because they form a semigroup under the 
action of the evolution operator for z > 0. [ 20 ] 

IV. REVERSED HARMONIC OSCILLATOR 


In the HNA, for a = 0, GV can be explicitly written: 
letting 


= 4- 

(10) 

we have. 


H = Pkq 

(11) 

with E^ = P /^0 5 and 


-2 2 2 


- p J X 

“ 2 2 ’ 

(12) 

with 7 ^ = P/^ 2 - For the exponential nonlocality, 
Kq = 1/2(7 and = I/a/Pct^. 

We let 0 = exp (^—ikqPz) 0, with 

we have 

o 

II 

(13) 


Hro is the Hamiltonian of the reversed harmonic oscil¬ 
lator (R-HO). [15, 25, 26] In the classical limit, the tra¬ 
jectories of the RHO obey x{t) = — 7 x(t), i.e., they cor¬ 
responds to a dissipative system, and were also found 
in the hydrodynamical approximation of DSW. [27] This 
shows that the nonlocal optical propagation in a defo- 
cusing medium is a physical realization of a quantum 
dissipative system. [28] GVs for the RHO are 

= (y - k = Etfn (14) 

with purely immaginary = ±i{n P |) with [29] 

ftix) =—^==Hn{\/Sjx)exp{Ttp^), (15) 

being Hn{x) the Hermite polynomials. For any n, two 
GV exist: displaying an exponential decay, and = 

(f“)* with exponential growth, and 

(fmifn) ~ (16) 

being 5mn the Kronecker symbol. is the time-reversed 

correspondent of f~. 

In Fig.la we show the profiles of the lowest order even 
GV; their modulus squared increases with x. We have 

oo 

Hno = E^n\fn){fi\ ( 17 ) 

n=0 

, and correspondingly 

oo 

4>^ix) = ( 18 ) 

n=0 

In the RHS theory the initial state '0((r,O) is prepared 
before the interaction, the state in the future is expressed 
by the that decays upon evolution. 

V. NONLOCAL SOLITONS AND DISPERSIVE 
SHOCK WAVES 

In the HNA, the fundamental soliton corresponds to 
the Gaussian shaped ground state of the HO.[ 22 ] The 
analytical prolongation to imaginary-valued energies of 
the HO are the GV. There is hence a precise link be¬ 
tween intrinsically irreversible states (the GV) and the 
propagation invariant self-trapped beam (the fundamen¬ 
tal soliton): the prolongation of the ground state of the 
standard HO (the fundamental soliton) to immaginary 
eigevalues, furnishes the fundamental GV fo (^)- 

There is also another notable connection with DSW. 
Not only GVs are found in the same regime of the non¬ 
local DSW; [30] but the phases of the GVs have the sig¬ 
natures that are commonly ascribed to DSW. In Fig.Ib 
we show the GV phase profile when increasing n, which 
progressively resemble DSW (e.g., compare with Fig. I of 
[30]). In these terms the singularity that arises upon the 
generation of shock waves, when increasing the power, 
can be interpreted as the excitation of higher order GVs. 
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VI. DECAY PROCESSES 

In the probabilistic interpretation of RHS-QM, the 
projection Eq.(18) over gives the probability 

Pn{z) of finding the system in a decaying GV. This gives 
the evolution of the irreversible part of the wavefunction, 
which can be also interpreted as the decaying process of 
elementary excitations with energy given by = KqP 
and quantized decay rates T^ = j{2n 1): 

N 

n 

The expansion in Eq.(19) is only meaningful for z > 0 
and 

Pn{z) = r„|(f+|^(a7,0))|^e"’"""= (20) 


gives the z—dependent weight of the n-order GV. If 
0) is a pure GV, 


II 

(21) 

we have 


Pn{z) = r„exp(-r„2;), 

(22) 

with normalization[26] 


POO 

/ Pn{z)dz = 1. 

Jo 

(23) 

The excited GVs depend on the initial profile. 
Gaussian beam. 

For a 

'0(x, 0) = i)F{x) = exp(—x^/2)/^^. 

(24) 

all the odd terms in Eq.(19) vanish due the x- 
and for the first two even GV we have 

-parity, 

^3/2 

Po{z)=‘^ /-- 2exp( 7z), 

Vl + 7 

(25) 

and 


r 7^/^ 

P 2 {z) = 5 -^ exp( 572;). 

VI+ 7^ 

(26) 

In Fig. Ic we show the weights Pn(0) for various GV 
versus the normalized power 7^ = ^P\ they tend to 
zero for P ^ 0 and grow with P. Their quantized decay 
rates scale with power as 7 oc a/P. 


We underline that GV are not decay processes depending 
on the coupling with the environment (i.e., extrinsic)^ 
but exponentially decaying states arising from the time- 
reversible Hamiltonian {intrinsic origin). 

To show the occurrence of the GVs in the original non¬ 
linear model, we solve Eq.(3) with '0(x,O) = ippix). At 
low power (not reported) no exponential decay is found. 
At high power (high wave amplitude), in correspondence 


of the DSW, the resulting dynamics clearly display expo¬ 
nential decays (Eig. 2a). The shape of the beam strongly 
resembles the excitation of the ground state GV, cor¬ 
responding to a central plateau; and lateral tails can 
be identified with higher order GVs (compare Eig. 2a 
with Eig. la). As the GVs decay exponentially and the 
power is conserved, the beam displays a self-similar evo¬ 
lution, with an exponential spreading following the clas¬ 
sical trajectories of the dissipative system. To provide 
quantitative evidence of , we project the wavefunc¬ 
tion 2 p{x^z) over and retrieve Pn{z) in Eq.(20), as 
shown in Eig.2b,c. Pn{z) decays with quantized rates 
Tn = (2n + 1)7. This trend has been verified for various 
7^ and initial conditions, and confirms the theoretical 
analysis. The evidence of the quantization of decay rates 
is the most direct signature of the GVs, as shown in Eig. 
2b,c. Deviations from the exact exponentials are due to 
the finite degree of nonlocality in Eq.(3). 


VII. GAMOW VECTORS AND THE DEGREE 
OF IRREVERSIBILITY 

DSW are formally reversible, but in real world exper¬ 
iments there is always an amount of loss. The interplay 
between linear losses and the excitation of the GV is sub¬ 
tle. We quantify in the following the way irreversibility, 
in the presence of linear losses, is enhanced by the exci¬ 
tation of nonlinear GV. In order to give a quantitative 
analysis, we define below the degree of irreversibility as 
the scalar between the forward and backward evolution 
in way such that when this product is zero the dynamics 
is totally irreversibile. We show below that the degree 
of irreversibility is enhanced by the nonlinearity with re¬ 
spect to linear loss, and this can quantitative calculated 
by retaining into account the excitation of GV, and using 
their exponential dynamics in the analysis. 

We start showing in Fig. 3a-d the numerical solutions 
of Eq.(3). We let '0(x,O) = i1^f{x) and retrieve '0(x,I/). 
We then use the conjugated propagated field as new ini¬ 
tial condition '0'(x,O) = '0(x,I/)*. As long as dynamics 
are reversible the propagated 'ipsix) = '0'(x, L) coincides 
with 'iprix) (Fig. 3a,b). Repeating the procedure for 
a = 0.3, reversibility does not occur (Fig. 3c,d). The 
discrepancy between and increases with P, for 
a > 0. We introduce the degree of reversibility by 

i?=KV'F|V'B)p. (27) 

For a fixed L, P = R{a, P), and in the absence of loss 
a = 0, the dynamics is exactly reversible 

P(0,P) = 1, (28) 

as show in the figure 3a,b. For a > 0, in the absence 
of nonlinear effects P = 0, the degree of reversibility is 
given by the linear loss 

P((a,0) = exp(—aP), 


(29) 
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increasing linear losses reduces reversibility. We calculate 
R{a, P) and, as shown Fig.Se, the breaking of reversible 
character appears more pronounced when increasing P. 
To discriminate extrinsic (only due to the loss a > 0) 
and intrinsic contribution to irreversibility (due to GV) 
we introduce the fraction of the intrinsic irreversibility 

Iint{oL^P) = 0 if the irreversibility is only due to linear 
losses; if P) > 0 nonlinearity enhances irreversibil¬ 

ity and lint quantifies this contribution. 

Figure 3f shows that Iint{oL^ P) grows with P. To pro¬ 
vide evidence of the role of GVs, we consider the evolu¬ 
tion according to Eq.(3) of the truncated (finite power 
realization) states, i.e., ' 0 (x,O) = f^(x) for \x\ < xq 
and zero elsewhere. In Fig.4a,b we show the dynamics 
of truncated decaying states n = 0 , 2 , strongly resem¬ 
bling that observed for a Gaussian initial condition at 
high power. Fig. 2 b, mediating the DSW generation. 

In Fig.4c,d we show the backward evolution of the ex¬ 
ponentially amplified truncated states f+, which are the 
time-reversed f“. Albeit no amplification is present, but 
also losses are included, these undergo a strong field en¬ 
hancement in propagation and collapse-like dynamics. 
This is due to the fact that the power conservation, for 
negligible loss, forces the exponential shrinking of the 
wavepacket spatial-width when the peak suffers the ex¬ 
ponential increase. The collapse of the time-reversed GV 
state is ultimately limited by the nonlocality. The dy¬ 
namics of time-reversed GV reproduces that during re¬ 
versibility breaking in Fig. 3d; i.e., the forward evolution 
generates the states that, upon time-reversal, becomes 
the exponentially amplified and collapsing = (f“)*. 

To explain in simple terms the interplay between the 
linear loss and the GV, we observe that after time re¬ 
versal the exponentially decaying dynamics of GV, with 
exponent — r^(P) — a, becomes an exponential amplifi¬ 
cation of states growing with r^(P') — a; on the con¬ 
trary, the background component (ftij is always attenu¬ 
ated. P' is in general different from P because of linear 
loss. The presence of GV radically alters the degree of 
reversibility because of their exponential dynamics: in a 
forward/backward round of evolution GV are nearly un¬ 
affected by linear loss, while the background component 
is attenuated. 

The previous arguments allow to derive a simple ex¬ 
pression for the degree of irreversibility: at high power 
the beam is mostly formed by GV; for the n—state we 
have 'iPf{x) = fn~(:^^)- After forward propagation, letting 
Leff = 2(1 — we have 

= fn“{a:)exp[-a^ - (31) 

Upon time reversal we have 

= fn~'~(x) exp[-a^ - TniP)^^], (32) 



FIG. 3. (a) forward progagation of 'iprix) for a = 0; (b) 

backward propagation corresponding to panel (a); (c) as in 
(a) for a = 0.3; (d) backward propagation corresponding to 
panel (b); (e) degree of reversibility versus P for various a; 
(f) fractional degree of intrinsic irreversibily; the dashed line 
is the estimated after Eq.(34). Parameters: L = 2, cr^ = 10 


which propagates again in the forward direction, at a 
power P' = Pexp(—(aP): 


by projecting over we have for n = 0 


I in 


1 — exp 



0 


-^e// 

2 


(34) 


which is found in good agreement with the numerical 
values in Fig.3e, also taking into account that we are 
only retaining the n = 0 term in the GV expansion. In 
the limit of small loss a, Eq.(34) reads as 


^/PaF y? / |n2|/o L\ for\ 

4 ^<7 4 V no LnlocLloss' ^ ’ 

being and Jq = Pmks the peak intensity Lz = 

LZd the propagation length in real world units. 

The use of GV radically simplifies the treatment of the 
highly nonlinear and nonlocal regimes. Eq.(35) shows 
that the strength of nonlinearity has a direct effect on 
irreversibility, due to states with intrinsic exponential de¬ 
cay. A small amount of extrinsic loss leads to a breaking 
of time-reversal that is amplified by nonlinearity. In these 
terms GV have a leading role, as demonstrated by this 
theoretical analysis that heavily relies on the use of GV 
exponential trend, which is shown to be in quantitative 
agreement with numerical simulations. 


VIII. CONCLUSIONS 

We theoretically predicted Gamow vectors in nonlin¬ 
ear waves and discuss their possible first observation in 
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FIG. 4. (a) forward propagation of the truncated GV 

(n = 0, xg = 40, P = 4, a = 0 , cr = 10); (b) as in (a) for 
n = 2 ; (c) backward propagation of truncated GG (n = 0 , 
XG = 40, P = 4, a = 0, cr = 10); (d) as in (c) for n — 2. 


laser beam propagation in thermal liquids. We discussed 
their main signature, i.e., the quantization of their decay 
rates, and show the strict link between Gamow vectors 
and nonlinear wave physics. 


In particular, we investigated the degree of irreversibil¬ 
ity of dispersive shock waves in highly nonlocal media. 
We show that the exponentially-decaying dynamics of 
Gamow vectors enhances the irreversibily, and makes it 
impossible reverting a dispersive shock wave even in the 
presence of a neglible and arbitrarily small loss. The ev¬ 
idence for this role of Gamow vectors is given by the fact 
that they allow the analytical calculation of the degree of 
irreversibility in quantitative agreement with numerical 
solutions. This is a subtle result of the interplay between 
the exponential decay of nonlinear origin and the cou¬ 
pling with environment. 

This analysis shows that a nonlocal nonlinear propa¬ 
gation may be used to simulate quantum systems and 
experimentally study the origin of irreversibility. Our re¬ 
sults are relevant for other problems in nonlinear physics 
such as rogue waves, beam collapse and filamentation, 
wave-breaking and supercontinuum generation, Bose- 
Einstein condensation and polaritons, and also provide 
support to the RHS formulation of Hamiltonian theories. 
The exponential dynamics of Gamow vectors may also 
potentially employed in applications as lasers and nano- 
focusing of light. 

We acknowledge support from the CINECA under 
the ISCRA initiative, the Progetto di Ateneo Sapienza 
Award: PhotoAnderson, and the ERG project VAN¬ 
GUARD, grant number 664782. 
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